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Abstract. The first-principles calculations have been performed to investigate the ground 
state properties of cubic scandium trifluoride (SCF3) perovskite. Using modified hybrid 
exchange-correlation functionals within the density functional theory (DFT) we have 
comprehensively compared the electronic properties of SCF3 obtained by means of the linear 
combination of atomic orbitals (LCAO) and projector augmented-waves (PAW) methods. 
Both methods allowed us to reproduce the lattice constant experimentally observed in cubic 
SCF3 at low temperatures and predict its electronic structure in good agreement with known 
experimental valence-band photoelectron and F Is X-ray absorption spectra. 



1. Introduction 

Scandium trifluoride (SCF3) is a promising material attracted much attention due to 
recently discovered strong negative thermal expansion (NTE) coefficient over a wide range of 
temperatures from 10 to 1100 K [1]. 

At atmospheric pressure SCF3 has cubic (space group Pm3m) Re03-type structure (Fig. [1]) 
down to at least 10 K pQ. However, opposite to the metallic Re03, which has very week NTE 
effect [21 [3] , scandium trifluoride is an insulator with band gap of more than 8 eV [4 J . A cubic- 
to-rhombohedral phase transition occurs in SCF3 at high pressure (P>0.5 GPa at T~300 K or 
P=0.1-0.2 GPa at 50 K) as determined from X-ray and neutron diffraction studies [HE]. Raman 
spectroscopy confirms that upon pressure increase cubic SCF3 undergo a phase transition to the 
rhombohedral (space group R3c) phase [6] . 

The electronic structure of SCF3 thin films has been investigated by means of resonant 
photoemission spectroscopy at the Sc 2p and F Is absorption edges in Ref. [4j. Two charge- 
transfer-type satellites, the first one at about 13 eV below the main peaks in the valence-band 
and the second one with an energy separation of about 9-10 eV, were observed and indicate the 
strong hybridization effect between Sc 3d and F 2p states [3J. 

The phonon properties of SCF3 have been studied using a combination of inelastic neutron 
scattering experiment with ab initio calculations of lattice dynamics [7J. It was shown that a 
description of NTE within the quasi-harmonic approximation is not reliable in the case of SCF3. 
The authors in [7] have demonstrated that the R4+ mode (the one with the lowest energy at 
i?-point of Brillouin zone (BZ)) has quartic potential and proposed the mechanism for NTE 
based on this anharmonicity. The transition between the ground state and the first excited 




Figure 1. Schematic view of cubic SCF3 
structure (space group Pm3m) built up of 
ScF6 regular octahedra joined by corners 
(small balls are fluorine atoms, large balls 
are scandium atoms). All Sc-F-Sc angles 
are equal to 180°. 



state (~19 meV) of the -R4+ mode has been associated with experimentally observed phonon 
peak (~25 meV), which stiffens upon increasing temperature [7]. 

In spite of mainly experimental efforts have been performed to understand the NTE in SCF3, 
its origin is still under debate. Therefore, comprehensive theoretical study is necessary to 
elucidate the nature of this peculiar property. In this study we report the results of ah initio 
electronic structure calculations for cubic SCF3 using the linear combination of atomic orbitals 
(LCAO) and projector augmented-waves (PAW) methods. 

2. Computational details 

In this study SCF3 is modeled by means of two different methods: (i) LCAO within the framework 
of hybrid density functional approach and (ii) PAW calculations using the generalized gradient 
approximation (GGA) density functional. 

To perform hybrid LCAO calculations, we used the periodic CRYSTAL09 code [8], which 
employs Gaussian-type functions centered on atomic nuclei as the basis set for expansion of 
the crystalline orbitals. The full-electron basis sets used in this study for F and Sc were taken 
from CRYSTAL basis set library [SJ. Threshold parameters for evaluation of different types of 
bielectronic integrals such as overlap and penetration tolerances for Coulomb integrals, overlap 
tolerance for exchange integrals and pseudo-overlap tolerances for exchange integral series [5J 
have been set to 10 -8 , 10~ 8 , 10 -8 , 10 -8 and 10 -16 , respectively. Calculations are considered as 
converged when the total energy obtained in the self-consistent field procedure differs by less 
than 10 -10 a.u. in the two successive cycles. 

LCAO calculations of SCF3 have been performed using a set of hybrid Hartree-Fock/Kohn- 
Sham (HF/KS) exchange-correlation functionals [9] 110 1 fTT |, [12 ] [13] combining exact HF nonlocal 
exchange and KS exchange operator within the GGA as implemented in the CRYSTAL09 code 
[HJ. We have constructed the following hybrid functionals: "BECKE & LYP", "BECKE & 
PWGGA" , "PBE & PBE" , "PBEsol & PBEsol" , where the first terms refer to the KS exchange 
and the second terms to the KS correlation parts of a hybrid functional, while the weight of the 
HF exchange has been varied for each functional separately between and 100%. 

DFT calculations using the PAW method as implemented in the VASP code [14] were 
performed for comparison. We have compared PW91 [13] and PBE [9] exchange-correlation 
potentials as well as large-core effective core potential (LC-ECP) with three 3d 1 4s 2 valence 
electrons and small-core effective core potential (SC-ECP) with eleven 3s 2 3p 6 3d x 4s 2 valence 
electrons for Sc atom. Fluorine atom for all calculations was described as ECP with 2s 2 2p 5 
valence shell. Hybrid HF-DFT calculations have been also performed with the weigth of exact 
HF exchange part equal to 25%. The optimal cut-off energy has been equal to 1000 eV. 

To provide the balanced summation over the direct and reciprocal lattices of SCF3 in both 
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Figure 2. Equilibrium lattice con- 
stant with respect to the weight of ex- 
act HF non-local exchange in hybrid 
density calculation by means of the 
LCAO method. 



LCAO and PAW methods, the reciprocal space integration has been performed by sampling 
the BZ with the 8x8x8 Pack-Monkhorst A;- mesh [15] that results in 35 eventually distributed 
/c-points in the irreducible BZ. 

3. Choice of exchange-correlation functional 

The influence of the exchange-correlation functional type and the weight (wHF) of exact HF non- 
local exchange in hybrid density calculations within the LCAO method on the lattice parameter 
in SCF3 is shown in Fig. [21 We found that pure GGA functionals overestimate the experimental 
value of lattice constant ao=4.026 A [lj, while stepwise increase of the u>HF for about 10% 
leads to proportional decrease of ao by approximately —0.01 A. The PBEsol functional with 
u>HF=18% reproduces perfectly the experimental lattice parameter and has been selected for 
further calculations. 

Values of the lattice constant ao obtained using pure GGA calculations by the VASP code 
[T3] are in the range of 4.02-4.03 A if LC-ECP is used for Sc, while a of 4.06-4.07 A has been 
obtained if SC-ECP has been adopted for Sc. Using hybrid HF-DFT calculations with the PBE0 
functional [H]), the value of ao is equal to ~4. 04-4. 05 A: it is larger when LC-ECP is used for 
Sc and is smaller when SC-ECP is used. 

To additionally justify the choice of hybrid functional, we have calculated the value of the 
band gap (E g ). Until now the band gap of SCF3 was not precisely measured experimentally, but it 
is estimated to be larger than 7-8 eV All hybrid functionals constructed for LCAO method 
(u>HF corresponds to ao=4.026 A, see Fig. [2]) show that SCF3 is an insulator with relatively 
large band gap E g > 8 eV. In the PAW GGA calculations the band gap is underestimated 
as it is predicted by the theory for pure GGA functionals: it is equal to 5.5-6.0 eV. Hybrid 
PAW calculations allowed us to obtain E g equal to 8-9 eV, which is closer to the experimental 
estimate. 

Thus, for electronic structure calculations we have finally chosen the "PBEsol & PBEsol" 
functional with iuHF=18%, giving E g =9.8 eV, for LCAO method and the PW91 functional 
with wHF=25%, giving E g =8.9 eV, for PAW method. 

4. Electronic properties 

The effective charges of scandium and fluorine atoms were calculated using Mulliken population 
analysis [8] in LCAO method and Bader topological analysis [16] in PAW calculations. The 
LCAO charges are equal to gs c =+2.28e and of=— 0.76e. In PAW method the atom charges 
were slightly different for two Sc pseudopotentials: (/sc=+2.70e and qp=— 0.90e for LC-ECP and 



Figure 3. Plots of the charge 
density difference in (100) plane 
between the full charge density 
of crystal and electronic density 
of non-interacting atoms (left 
panel) and ions (right panel) ob- 
tained using "PBEsol & PBEsol" 
functional with wHF = 18%. 
Solid (red) and dashed (blue) iso- 
lines correspond to positive (ex- 
cess) and negative (deficiency) 
electron density difference, re- 
spectively. 

Figure 4. A comparison of the 
F p density of states (DOS) cal- 
culated by means of the LCAO 
(dashed line) and PAW (solid 
line) methods with valence-band 
photoelectron [17] and F Is X- 
ray absorption spectra [I]. The 
amplitude of the DOS in the con- 
duction band was multiplied by 
a factor of 5. Energy scale is 
given with respect to the calcu- 
lated top of the valence band. 
Experimental spectra has been 
aligned to match best the calcu- 
lations. 

c/Sc=+2.17e and qF=— 0.73e for SC-ECP. These results demonstrate a considerable deviation of 
effective atomic charges from formal charges caused by partly covalent nature of Sc-F bonds. 

Fig-E]shows electron density maps in (100) plane of SCF3 obtained using LCAO calculation as 
a difference between the full charge density of crystal and electronic density of non-interacting 
atoms or ions. The anisotropy of electron density induced by bond covalent contribution is 
clearly visible in both graphs: the electron density is larger along the crystallographic axes. 

In Fig. U] we compare the partial F p density of states (DOS) calculated by both LCAO and 
PAW methods with the experimental valence-band (VB) photoelectron spectrum taken from 
Ref. [T7] and the F Is X-ray absorption spectrum taken from Ref. [I]. The VB photoelectron 
spectrum of SCF3 has a complex shape composed of the main peak and the shoulder located at 
the left side, which originate due to the valence band splitting upon hybridization of the F 2p 
and Sc 3d states [17] ■ Our theoretical calculations confirm such interpretation and reproduce 
well the width and asymmetric shape of the valence band. The F Is absorption spectrum [3] 
consists of three well visible peaks at 9.3, 11.4 and 18.2 eV due to transitions from the F Is core 
state to the F np states in the conduction band and above. The first two peaks correspond to the 
F 2p states in the conduction band, which are hybridized with the Sc 3d states. The two peaks 
originate from the t2 g -e g splitting of the Sc 3d states by the octahedral crystal field. The third 
peak at 18.2 eV corresponds to the interaction of the F 2p and Sc 4p states. Such interpretation 




is fully supported by the calculated F p DOS: in fact, good agreement with the F Is absorption 
spectrum indicates that a relaxation of the electron states due to the F Is core- hole is weak in 
ScF 3 . 

5. Conclusions 

The lattice constant, band gap and electronic properties of cubic SCF3 were calculated using 
LCAO and PAW DFT computational methods. It was shown that a variation of exact HF 
exchange part allowed us to reproduce experimental value of the lattice constant ao=4.026 A 

m- 

The band gap obtained in hybrid HF-DFT calculations is equal to 8-10 eV and agree well 
with the expected value [1J. A deviation of atom charges from formal ones and an anisotropy 
of electronic density point to the covalent nature of Sc-F bonds in SCF3. Calculated partial F 
p DOS is in a good agreement with valence-band photoelectron [17] and F Is X-ray absorption 
[I] spectra. 

Good agreement between the theory and available experimental data makes promising the 
application of the present approach to the investigation of phonon properties in SCF3 in the 
future work. 
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